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ABSTRACT 

We present kinematic measurements of a large sample of galaxies from the Team Keck Redshift Sur- 
vey in the GOODS-N field. We measure line-of-sight velocity dispersions from linewidths of integrated 
emission for 1089 galaxies with median redshift 0.637, and spatially resolved kinematics for a subsam- 
ple of 380 galaxies. This is the largest sample of galaxies to z ~ 1 with kinematics to date, and allows 
us to measure kinematic properties without morphological pre-selection. Emission linewidths provide 
a dynamical measurement for the bulk of the blue galaxy population. To fit the spatially resolved 
kinematics, we construct models with both line-of-sight rotation amplitude and velocity dispersion 
as fit parameters. Integrated linewidth correlates well with a combination of the spatially-resolved 
velocity gradient and dispersion, and is a robust measure of galaxy kinematics. The spatial extents of 
emission and continuum arc similar and there is no evidence that linewidths are affected by nuclear or 
clumpy emission. The measured rotation gradient is a strong function of slit position angle alignment 
with galaxy major axis, but integrated linewidth is largely independent of slit alignment. Even in 
a subsample of galaxies with well- aligned slits, there are galaxies whose kinematics are dominated 
by dispersion (V/a < 1) rather than rotation. These are probably objects with disordered velocity 
fields, not dynamically hot stellar systems. About 35% of the spatially resolved sample are dispersion 
dominated; galaxies that are both dispersion dominated and bright exist at high redshift but appear 
rare at low redshift. This kinematic morphology may yield a probe of the evolutionary state of these 
objects. Kinematic morphology is linked to photometric morphology in HST/ACS images: dispersion 
dominated galaxies include a higher fraction of irregulars and chain galaxies, while rotation dominated 
galaxies are mostly disks and irregulars. Only one-third of chain/hyphen galaxies are dominated by 
rotation; high redshift elongated objects cannot be assumed to be inclined disks. In a companion 
paper, we use the linewidths and rotation to measure evolution in the Tully-Fishcr relation. 
Subject headings: galaxies: distances and redshifts — galaxies: evolution — galaxies: fundamental 
parameters — galaxies: high-redshift — galaxies: structure — surveys 



1. INTRODUCTION 

Galaxy redshift surveys to z ~ 1 allow the measure- 
ment of evolution in properties of the galaxy population 
at lookback times of half the present age of the Uni- 
verse. To date, a number of surveys have studied galaxy 
properties such as luminosities, colors, radii, morphology 
and environment. The internal kinematics of galaxies arc 
equally interesting, as a physical probe of individual ob- 
jects and to allow the measurement of scaling relations of 
the galaxy population. However, measuring kinematics 
requires moderate to high-resolution spectroscopy, which 
has until recently been difficult to obtain for large num- 
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bers of faint galaxies. 

This paper measures kinematics of galaxies from emis- 
sion lines, using Keck/DEIMOS spectra from the Team 
Keck Redshift Survey (TKRS) in the GOODS-N (Great 
Observatories Origins Deep Survey) field (Wirth et al. 
2004; Giavalisco et al. 2004). The TKRS provides spec- 
tra of 1437 galaxies drawn from a magnitude limited 
sample to Rab < 24.4, at resolution R ~ 2100. We 
select 1191 galaxies with good emission line detections 
and measure 1089 emission line velocity dispersions from 
integrated emission, with median < z >= 0.637. A sub- 
sample of 464 are selected for modeling spatially resolved 
rotation and dispersion profiles, of which 380 yield good- 
quality measurements. We use these spatially resolved 
measures to probe the kinematic properties of high red- 
shift galaxies and to test the integrated linewidths. In 
Paper II (Weiner et al. 2006), we use the linewidths 
and rotation velocities to measure evolution in the Tully- 
Fisher relation. 

A number of previous works have observed galaxy in- 
ternal kinematics from 0.1 < z < 1.0, for the purpose of 
measuring Tully-Fisher relations. The pioneering studies 
of Vogt et al. (1996, 1997) modeled rotation curves for 17 
galaxies of disky morphology with median < z >= 0.47 
by combining Keck/LRIS slitlet spectra with structural 
information from HST photometry. Several subsequent 
studies of rotation curves with similar modeling proce- 
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dures have contained 20-100 galaxies with median red- 
shifts - 0.4 - 0.5 (e.g. Simard & Pritchet 1998; Vogt 
2000; Ziegler et al. 2002; Milvang- Jensen et al. 2003; 
Bohm et al. 2004; Bamford et al. 2005,2006; Conselice 
et al. 2005; Nakamura et al. 2006; Metevier et al. 2006). 
Generally these samples have been morphologically se- 
lected to be inclined disky objects, of which a majority 
show measurable rotation. 

Here we study a sample with emission line kinemat- 
ics that is essentially selected only on magnitude and 
emission line strength. We measure one-dimensional in- 
tegrated velocity dispersion (linewidth), and for a sub- 
sample, spatially resolved rotation profiles; we discuss 
the properties of these velocity measures in Sections 0] 
and |5l We show that the integrated linewidth is a fairly 
robust measure of the characteristic velocity of a galaxy, 
expanding the scope of the study beyond galaxies that 
are selected to be orderly rotating disks. Indeed, we find 
from the spatially resolved rotation profiles that a signif- 
icant number of high-redshift galaxies show kinematics 
that do not appear to be orderly rotation. 

We adopt a LCDM cosmology with h — 0.7, f2jvf = 0.3, 
and Q\ = 0.7. Magnitudes quoted in this paper are in 
the AB system unless explicitly indicated as Vega. Sec- 
tion |2 discusses the sample and measurement methods, 
and section presents the properties of the sample and 
completeness. Sections 0] and [S] study the simulated and 
empirical properties of the kinematic measurements: in- 
tegrated linewidths and spatially resolved velocities. In 
Paper II we use the kinematics to measure the evolution 
in the Tully-Fisher relation with redshift. 

2. SAMPLE, VELOCITY MEASUREMENTS, AND 
PHOTOMETRY 

2.1. Properties of the spectroscopic sample 

The parent sample for our kinematic measurements is 
the Team Keck Redshift Survey (TKRS) in the GOODS- 
N field. The selection and observations are described 
at length by Wirth et al. (2004) and we present a 
few relevant details here. Data and images of TKRS 
galaxies may be retrieved from the TKRS website at 
http : //www2 .keck. hawaii . edu/science/tksurvey/ 

The TKRS sample is magnitude selected to Rab = 
24.4. No color or morphological selection was applied, 
but there is a mild surface brightness selection (see 
Willmer et al. 2006), which may eliminate a few faint 
high-z reddened edge-on galaxies. The sampling com- 
pleteness (ratio of galaxies with spectra to galaxies meet- 
ing selection criteria) is ~ 75% and the redshift success 
rate is ~ 70% (> 80% for R AB < 23). Most galaxies for 
which spectra were taken but redshifts could not be de- 
termined are faint, blue, and probably at z > 1.5, where 
[O n] A 3727 disappears into the night-sky OH line forest. 

TKRS observations used the Keck II telescope and 
DEIMOS spectrograph with a 600 lines/mm grating, 
yielding a sampling of 0.119"/pixel and 0.648 A/pixel. 
Typical wavelength coverage was 4600-9800 A. 18 slit- 
masks were observed, each for 60 minutes total exposure 
time. The slit width was 1.0". In the mask designs, slit 
position angles (PAs) could be tilted to follow the ma- 
jor axis of elongated objects, by up to 30 degrees away 
from the nominal perpendicular-to-dispersion direction. 
Galaxy PAs were measured from ground-based i?-band 
imaging before mask design. When slits were tilted, the 



slit width was adjusted to keep the width along the dis- 
persion direction at 1.0", keeping the spectral resolution 
constant regardless of slit PA. 

The TKRS data were reduced and redshifts measured 
with the DEIMOS pipeline software developed by the 
DEEP project, described in Davis et al. (2003) and 
Wirth et al. (2004); TKRS data are very similar to 
DEEP2 data, save that DEEP2 uses a 1200 lines/mm 
grating. Candidate redshifts were measured automat- 
ically and verified by visual inspection by members of 
Team Keck. 

End data products include 2-d spectra for each slitlet 
and both optimal and boxcar extracted 1-d spectra for 
each object. We used the 2-d spectra to measure rota- 
tion curves and the 1-d spectra to measure linewidths 
of integrated emission. In this paper we always use the 
boxcar extraction, never the optimal: the optimal extrac- 
tion weights different regions of the galaxy light profile 
unequally, which makes the physical meaning of the op- 
timally extracted spectrum somewhat obscure. The box- 
car extraction window diameter is set by the reduction 
software to 1.5 times the FWHM of the object measured 
in the 2-d spectrum. For nearly all galaxies, larger ex- 
traction windows do not make a significant difference to 
the measured linewidths, because the emission intensity 
is centrally peaked (see also Section l4~2"|) . 

2.2. Emission line measurement in the 1-d spectra 

Our goal in fitting linewidths to integrated emission is 
to obtain a kinematic measurement for as many galaxies 
as possible, without applying cuts to the sample. Several 
previous studies have used linewidths of integrated emis- 
sion to measure kinematics of small samples of galaxies 
at intermediate redshift (Rix et al. 1997; Mallen-Ornelas 
et al. 1999). 

2.2.1. Line fitting 

We developed an automated program to fit emission 
lines in all TKRS galaxies with secure redshifts (quality 
code 3 or 4 indicating two spectral features, Wirth et al. 
2004). The LINEFIT program takes a galaxy redshift 
and list of common emission lines; at the predicted loca- 
tion of each line, it takes the wavelength, flux, and error 
data from a 40 A window and fits a gaussian profile, us- 
ing a Levenberg-Marquardt non-linear least squares x 2 
minimization (Press et al. 1992). 

The fit has four free parameters: continuum level, line 
intensity, velocity, and velocity dispersion. For the [O II] 
AA 3726.0, 3728.8 doublet, the program can fit a doublet 
with the wavelength ratio fixed and the intensity ratio 
fixed or free; for the TKRS data, we fixed the doublet 
intensity ratio at 1.4 (the mean ratio of red/blue compo- 
nents in high-S/N data measured with the 1200 lines/mm 
grating). The least squares fitter yields best-fit values 
and error estimates for all parameters. Fitting in the 
40 A window does not provide an adequate measure of 
the continuum in low-S/N spectra and so LINEFIT also 
measures a robust continuum by taking the biweight of 
data in two 80 A windows on either side of the emission 
line. 

2.2.2. Instrumental resolution 

This sample relies heavily on the measurement of ve- 
locity dispersion, which in turn can be strongly affected 
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TABLE 1 

Dispersion and velocity-related quantities 



Quantity units description 



by the instrumental spectral resolution (Ji ns t- Conven- 
tionally, the restframe intrinsic line-of-sight velocity dis- 
persion aid of a line is given by 

a ld = T J°l bs - °inst> CO 

where aid is in km s , A fc s is the observed wavelength 
in A, and a b s and ai ns t are the measured line disper- 
sion and instrumental resolution, both in A. As a rule 
of thumb, measurements for aid < c&inst/^obs are not 
very reliable because small errors in the observed width 
have a large effect on the inferred dispersion. We use 
aid to denote the restframe velocity dispersion derived 
from the 1-dimensional extracted spectrum, which is in- 
tegrated over the extraction window. Throughout these 
papers, observed quantities are given in A and restframe 
velocity quantities in km s _1 . This paper refers to sev- 
eral dispersion and velocity quantities, in observed and 
restframe, and 1-d and 2-d spectra; these are summarized 
in Table □ 

With DEIMOS and 1.0" slits, the profile of night sky 
lines or calibration arcs is somewhat flat-topped, less 
peaked than a gaussian. This flat-topping is pronounced 
for night sky lines in the 600 lines/mm grating data of 
the TKRS, but is relatively small for DEIMOS spectra 
taken with the 1200 lines/mm grating, such as those in 
the DEEP2 survey (Davis et al. 2003). Fitting gaussians 
to many line profiles of isolated night-sky lines in many 
different slits and masks yields mean instrumental gaus- 
sian sigmas ai nst ^k y — 1-4, 0.56 A for TKRS and DEEP2 
respectively. 

Flat-topping is less noticeable in the actual emission 
lines of observed galaxies in TKRS, for two reasons. 
First, galaxies' velocity widths broaden the profile so that 
flat-topping is smeared out. Second, TKRS galaxies are 
relatively small and their light distributions peak near 
the slit center, rather than uniformly filling the slit as 
night sky lines do. This slit underfilling is potentially a 
significant problem for any study of velocity dispersions 
of small objects. However, seeing and the relatively slow 
change of angular diameter with redshift beyond z ~ 0.5 
mitigate the effect. The effect is also smaller when the 
instrumental resolution is moderately high, as in TKRS 
and even more so for DEEP2. 

Modeling slit underfilling for the measured galaxy sizes 



and comparing objects with both 600 and 1200 lines/mm 
observations shows that flat-topping and its opposing ef- 
fects nearly cancel out. We matched observed linewidths 
for 70 galaxies on a DEEP2 slitmask that was observed 
with both 600 and 1200 lines/mm gratings, and found 
good agreement for 600 lines/mm ai nst = 1.4 A. 

The instrumental resolution as measured from night 
sky lines varies by ~ 5% peak-peak over the DEIMOS 
field: lines from slitlets near the ends of the mask are 
broader than lines from slitlets near the center. The ef- 
fect is only visible statistically, since its size is compara- 
ble to the scatter in individual sky line measurements and 
the scatter between different masks (focus variations). 
We have adopted a single value for the instrumental res- 
olution and neglected spatial and temporal variations. 
The variations' effect on the derived dispersions is small: 
for a small galaxy with aid = 50 km s , the peak error 
induced is ~ 2 km s , and the error declines rapidly for 
larger a xd - 

2.2.3. Velocity dispersion sample 

There are 1437 galaxies with redshifts, TKRS spectra, 
and magnitudes. We fit emission lines in each of these 
objects. For velocity dispersion purposes we then reject 
all lines that do not have a 4 sigma intensity detection. 
The 4 sigma line catalog contains 2595 lines over 1191 
galaxies. 

Galaxies can have several lines, and the fit parame- 
ters are independent. To obtain one estimate of velocity 
dispersion for each object, we take the weighted mean 
of the measurements of squared intrinsic velocity disper- 
sion, aid 2 ] using the square properly accounts for lines 
that are narrower than the nominal instrumental resolu- 
tion, which have aid 2 < 0. We use a weighted rms of 
the dispersions as the error estimate. We exclude a tail 
of 85 galaxies that have both error(log aid) > 0.25 and 
erior(aid) > 30 km s _1 to reject low-quality fits. 

Some line fits have an observed width a b s that is 
smaller than the nominal instrumental width. Although 
this is formally physically impossible, it is expected in 
the presence of noisy data, slit underfilling, and varia- 
tions in the instrumental resolution. After combining 
all lines, 196 of 1089 galaxies are "kinematically unre- 
solved," with widths close to or less than instrumental, 
so that the formal velocity dispersion is undefined or has 
a large error in log aid- When we restrict to Mb < —18, 
the cutoff for our Tully-Fisher fits in Paper II, 104 of 913 
are kinematically unresolved in log aid, under the crite- 
ria error(log aid) > 0.25, error(crid) < 30, and aid < 25 
km s _1 . Eliminating these galaxies from the sample pref- 
erentially rejects low-velocity galaxies and leads to a bias, 
so for plotting and fitting purposes we assign them a low 
value, log aid = 1.4 ± 0.2. The results of fitting do not 
depend strongly on the exact value assigned. Note that 
a few other galaxies have log aid < 1-4 yet low formal 
error on log aid, usually because they have very strong 
and well-measured emission lines. In Section 3 of Paper 
II we outline a fitting method which treats the unre- 
solved galaxies more robustly by fitting the ensemble of 
observed width a a b s before the instrumental resolution is 
subtracted. 

Some galaxies with very bright emission lines can have 
formally very small errors. It is hubris to take these 
errors literally, since we can hardly expect to measure 
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velocity dispersions of distant galaxies to < 5%. To be 
realistic and to prevent small errors from dominating fits, 
we add 0.03 in quadrature to the error on log aid, to cut 
off the low tail of the error distribution. The results are 
not sensitive to the exact value used; the median error 
for the sample of 913 brighter than Mb = —18 is 0.084 
in log aid- As discussed in Paper II, the intrinsic scatter 
in log aid dominates over the errors on individual points. 

2.3. Rotation curve measurements 

The DEIMOS slitlet spectra preserve spatial informa- 
tion and it is common to see emission lines with a veloc- 
ity gradient. Spatially extended emission can be used to 
measure a rotation curve. However, the spatial extent of 
emission is only a few times the seeing, so "beam smear- 
ing" is large and the effect of seeing must be modeled. 
The rotation curve is also affected by the inclination, slit 
width, and slit alignment with respect to the major axis 
of the galaxy. We measured rotation curves for a sub- 
sample of objects in order to study the relation between 
rotation velocity and integrated lincwidth. 

2.3.1. Previous high-redshift rotation curve modeling 

Previous efforts to measure rotation curves at interme- 
diate to high redshifts have approached these problems 
by modeling the intrinsic emission intensity and velocity 
distribution, subjecting them to seeing and slit effects, 
and fitting models to the data to obtain a measure of 
the rotation velocity. Vogt et al. (1996, 1997) pioneered 
work in this field, using HST/WFPC2 imaging to model 
the galaxy light distribution and parameters such as in- 
clination and position angle, and fitting to Keck/LPJS 
slit spectra, varying the model rotation velocity. Subse- 
quent works in this field have followed similar modeling 
programs, applying seeing and slit effects to a 2-d model 
of the light and velocity fields (e.g. Simard & Pritchet 
1998; Ziegler et al. 2002; Milvang- Jensen et al. 2003; 
Bohm et al. 2004; Bamford et al. 2005, 2006; Conselice 
et al. 2005; Nakamura et al. 2006; Metevier et al. 2006). 

These programs have generally focused on relatively 
small numbers of galaxies (~ 10 — 100), in some cases 
selected on morphology or spectra to be fairly normal ro- 
tating disks. They also require high-resolution imaging, 
preferably from HST, to obtain structural parameters. 
Most galaxies in the TKRS have deep HST/ ACS imag- 
ing from the GOODS survey (Giavalisco et al. 2004), 
and full modeling of the rotation curves is in progress. 
However, such modeling requires fitting of structural pa- 
rameters to ACS multidrizzled data and is beyond the 
scope of this paper. 

2.3.2. Rotation curve modeling from spectra alone 

In the spirit of obtaining kinematic measurements for 
as many galaxies as possible, we developed a simplified 
method for fitting seeing-compensated rotation curves 
using only the information contained in the 2-d spec- 
tra. This program, ROTCURVE, works by construct- 
ing models of the unblurred spatially resolved emission 
intensity, velocity and dispersion, blurring them along 
the slit to model seeing, and fitting to the data. Be- 
cause it only has one dimension of spatial information, 
the rotation model and the seeing convolution are 1-d 
functions along the slit, and there is no compensation 
for disk inclination or slit position angle. 



We chose a subsample of galaxies for rotation curve 
fitting using cuts in intensity of the strongest emission 
line and in spatial extent along the slit. The strongest 
emission line was required to have integrated intensity 
> 3000 e~ /pixel A in the summed 1-d spectrum. A 
gaussian was fit to the light profile along the slit and we 
required a Hght > 0.4", equivalent to FWHM > 0.94". 
464 galaxies were selected for ROTCURVE fitting. Of 
these 445 also had ACS imaging and ellipticity/position 
angle data, and 380 of the rotation curve fits were judged 
to be good by visual inspection. 

ROTCURVE fits to a single emission line for each 
object. It first fits an emission line to data in each row 
of pixels, using a similar algorithm to LINEFIT, to ob- 
tain profiles of velocity and dispersion along the slit, and 
rejects rows with discrepant values using automatic crite- 
ria, testing for e.g. large row-to-row jumps. It measures 
the light distribution along the slit of the continuum plus 
emission, fits a gaussian to this profile, and subtracts the 
assumed seeing (here 0.7") in quadrature to determine 
the intrinsic, unblurred intensity profile, modeled as a 
gaussian G{x) with dispersion r/. A minimum rj = 0.2" 
is imposed to keep intrinsic profiles of small objects from 
becoming pointlike. 

About half the TKRS masks have seeing — 0.7" and 
the remainder have seeing up to 1.0", determined from 
the spectra of identified stars and the distribution of 
fitted widths of galaxies. Assuming a smaller seeing 
than the actual leads to a small underestimate of V ro t- 
Assuming a larger seeing than the actual causes the 
ROTCURVE modeling to fail on small objects, because 
it derives a too-small intrinsic spatial extent 77, and the 
rotation velocity is poorly constrained due to the as- 
sumed large seeing blur. 

ROTCURVE then constructs models of the intrinsic 
position- velocity distribution along the slit. The veloc- 
ity model is an arctangent rotation curve centered at 
the peak of continuum+emission light, with a spatially- 
constant velocity dispersion, so that the intensity in 
position-velocity space before blurring by seeing and in- 
strumental resolution is: 



I{x,v) = G(x)exp(- 



v{x)f 



G(x) 



2nrj 



exp(- 



2a 2d 2 

(x - x Q y 
2rj 



V(x) = V rot — arctan(x/r t) ). 

7T 



(2) 



(3) 



(4) 



G{x) is the intrinsic gaussian light distribution along 
the slit. V(x) is the rotation curve with asymptotic cir- 
cular velocity V ro t and knee radius r v , and aid is the 
velocity dispersion, assumed constant along the slit. For 
the case of the [O 11] 3727 doublet, we use a double gaus- 
sian for the distribution in velocity, with intensity ratio 
1.4. 

The unblurred I(x, v) is then blurred for seeing with 
a 1-d gaussian in the spatial direction. The seeing has 
the effect of mixing gas at different velocities together, 
smoothing out velocity gradients and increasing the ob- 
served velocity dispersion at gradients (see Section l4~^l for 
an illustration of this effect). ROTCURVE then takes 
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moments of the blurred I i ur (x,v) in the velocity direc- 
tion to produce a model rotation curve and dispersion 
profile, and computes \ 2 of the data-model. It also fits 
for an offset in velocity since the published redshifts can 
vary from the systemic velocity, especially when there 
are large velocity gradients. 

The model parameters that can be varied are V ro t, (T2d, 
and r v . For each galaxy, we construct a grid of models 
over V ro t and <i2d- Because the models are fairly sim- 
ple and involve a series of 1-d convolutions, a brute-force 
minimization is practical; the minimum in \ 2 on the grid 
is always unambiguous and well localized. The best-fit 
values of 02d are near-quantized in units of the grid spac- 
ing (here 5 km s" 1 ) due to the minimum-finding tech- 
nique we used. 

The rotation scale length r v , describing how fast the 
rotation curve rises, is not strongly constrained by the 
ROTCURVE approach due to the seeing blur and the 
lack of 2-d modeling and position angle information. We 
chose a fixed r v = 0.2" for all galaxies, which can fit 
most rotation profiles; 0.1" and 0.3" fit worse for most 
galaxies. The knee radius and circular velocity are mod- 
erately covariant because the rotation curve turnover is 
not well resolved; changing r v from 0.2" to 0.1"or 0.3 
" changes V ro t by about 0.1 dex in the mean; larger r v 
causes larger V ro t- The effect on V ro t caused by assuming 
0.5" or 0.9" seeing is about 0.05 dex. The combination 
of V ro t and <T2d in quadrature, discussed in Section I5~2l 
below, changes by only half as much as V ro t- 

Figure ^ shows ACS /-band images of four example 
galaxies at z ~ 1 from the TKRS. Figurc[2shows postage 
stamp images of their emission line spectra, and Fig- 
ure shows their rotation and dispersion profiles fit by 
ROTCURVE. In Figure the points are velocity and 
dispersion from the spectra; the dashed and solid lines 
show the best-fit model before and after applying the 
seeing blur. These galaxies are relatively clean examples 
of a dichotomy in kinematic profiles that is common in 
the TKRS data. 

TKRS 2021 and 3168 show the familiar shape of a 
rotation curve. Beam smearing by the seeing not only 
smooths out the gradient in velocity, but induces a peak 
in the observed dispersion; the intrinsic dispersion is sig- 
nificantly less than the observed. TKRS 5627 and 10023 
are less familiar. They have quite small rotation, but 
require a velocity dispersion component to fit the data, 
with <j 2 d > V rot . Galaxies like TKRS 5627 and 10023 are 
common in our sample and are the reason we formulated 
ROTCURVE to include the velocity dispersion a 2 d as a 
second free parameter in the modeling. Note that these 
galaxies are relatively small but not extreme starbursts; 
they are not compact, nor very faint. We discuss rotation 
and dispersion dominated galaxies further in Sections l5.ll 
and below. 

We rarely see a rotation curve in which the observed 
points clearly roll over onto the flat part of the rotation 
curve, especially at z > 0.4. This is partly due to the 
limited extent of the data, but frequently a more impor- 
tant effect is the seeing blur, which makes the observed 
rotation curve rise shallower than the true rise. The ef- 
fect seen in TKRS 3168 in Figure |21 is common: the de- 
blurred model rotation curve has neared its asymptotic 
value while the observed data are still rising. Despite 
the lack of a flat-part in the observed data, the fitted 



values faithfully represent the observed velocity spread. 
Objects in which the data covered only a small spatial 
extent and the fitted velocity extrapolated beyond the 
data were rejected as erroneous fits. 

Because the rotation curve data points are highly cor- 
related, x 2 does not yield errors on V ro t and U2d directly. 
The fitted velocity V rot is a good representation of the 
gradient seen in the rotation curve, but subtler details 
such as the true errors on V ro t and its dependence on 
the spatial extent of the data will require further Monte 
Carlo simulation. Here and in Paper II, measuring Tully- 
Fisher evolution with V ro t and 02d is secondary to using 
them to probe the nature of galaxy kinematics, and as a 
test of the meaning of integrated linewidth a\d- 

As a check on the fits and models, ROTCURVE uses 
the observed row-by-row-fit rotation velocity and disper- 
sion profile to reconstruct the 2-d data, and subtracts it 
from the original data to produce a residual map. It also 
produces reconstructed 2-d data from the best-fit blurred 
model (convolving back to the instrumental resolution). 
The reconstruction from the row-by-row fits is generally 
very good. The reconstruction from the model is also 
generally good, although it leaves residuals in cases such 
as an asymmetry in emission intensity, since the model 
is symmetric by construction. In principle, the models 
could be fit directly to the 2-d data without ever reduc- 
ing to a 1-d rotation curve, but the rotation curve is a 
useful tool for evaluating the data and model. 

The ROTCURVE fitting procedure runs automati- 
cally. Afterwards we inspected the rotation and disper- 
sion profiles to reject bad or erroneous fits. 380 of the 445 
galaxies passed this inspection. The dominant causes of 
failures are: discrepant points in the rotation curve which 
evade ROTCURVE's automatic rejection; asymmetry 
which offsets the photometric center from the kinematic 
center; or small spatial extent of the rotation curve: we 
rejected most objects for which the usable rotation curve 
data spanned less than 1". 

ROTCURVE departs from previous high-redshift ro- 
tation curve modeling by making the velocity dispersion 
(J2d a free parameter; to our knowledge, other models 
have all fixed the dispersion to some low value. Fixing 
the dispersion is probably acceptable for samples of only 
morphologically-normal inclined disk galaxies observed 
with well-aligned slits, where rotation is the chief kine- 
matic support and is unambiguously detected. However, 
as seen for TKRS 5627 and 10023 and shown in Section 
15.11 there are galaxies with low rotation and large disper- 
sion in our broader sample, and for these, fitting without 
a dispersion term will yield erroneous rotation velocities. 

The dispersion term <j2d in ROTCURVE does not 
have to correspond to a literal gas velocity dispersion 
similar to the velocity dispersion of stars in an spheroidal 
galaxy. Because we are measuring nebular emission lines 
that occur in gas at T ~ 10 4 K, the gas in any individual 
H II region will have a dispersion of only 8-10 km s _1 . 
Taken literally, a dispersion &2d ~ 30 km s _1 implies a 
gas temperature of T ~ 10 5 K, which is an unfavorable 
place on the ISM cooling curve. It is more likely that the 
dispersion comes from relative velocities of discrete H II 
regions. 

Values of <j2d ^ 20 km s _1 can represent an effective 
dispersion, caused by the blurring of velocity gradients on 
scales at or below the seeing limit. Such a seeing-induced 
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Fig. 1. — ACS /-band images for the example galaxies whose rotation and dispersion are shown in Figure 131 TKRS 2021, 3168, 5627, 
and 10023. The images are 3" on a side; north is up and west is to the right. The panels are labeled with redshift (in the image), TKRS 
ID, restframe Mg , U - B, and linewidth log a ld . The DEIMOS slit PAs are 51°, 14°, 37°, and 39° respectively, north through east. The 
galaxies were classified as edge-on, LSB/irregular, chain, and hyphen/chain morphology, respectively. 




7460 7470 
wavelength 



i- 




8470 
0010023 



8480 8490 
wavelength 




Fig. 2. — Postage stamps of the DEIMOS spectra for the example galaxies whose rotation and dispersion are shown in Figure|3 TKRS 
2021, 3168, 5627, and 10023. The spectra have been sky-subtracted, and rectified for display purposes so that columns are constant 
wavelength. The galaxy continua, at low S/N/pixel, are not visible in these postage stamps. Each image shows the emission line used 
for ROTCURVE fitting, the [O n] 3727 doublet for 2021, 3168, and 10023, and [O III] 5007 for 5627. TKRS 2021 and 3168 show tilted 
emission lines, the signature of rotation curves, but TKRS 5627 and 10023 show very little velocity gradient. 



dispersion o~ia, can appear for slits which are aligned close 
to the minor axis of the galaxy and miss the resolved ve- 
locity gradient, as shown in Section FOl But for TKRS 
5627 and 10023, and similar galaxies in our sample, the 
slit was well aligned with the photometric major axis, 
yet there is almost no ordered velocity gradient. The ef- 
fective dispersion aid must be explained by non-ordered 
or disturbed gas kinematics on angular scales well below 
the seeing blur; these motions may not even have a pre- 
ferred plane. The geometry of the H II regions, and why 



they have not dissipated into a rotating disk, cannot be 
measured with these seeing-limited observations. 

2.4. Photometry 

Nearly all TKRS galaxies have BViz photometric 
measurements from the GOODS ACS catalog (Gi- 
avalisco et al. 2004). Ground-based photometry in 
UBVRIz'{HK') for a large field including GOODS-N 
is available from Capak et al. (2004). We used the 2.8" 
diameter aperture magnitudes from the GOODS cata- 
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(c) (d) 

Fig. 3. — Rotation and dispersion profiles for four example galaxies with well aligned slits. TKRS 202f and 3f68 at z = 1 .00 and 1.28 
are rotation dominated; TKRS 5627 and 10023 at z = 0.85 and 1.02 are dispersion dominated. The slit PA offsets are 6°, 9°, 3°, and 
10°. The points are the observed rotation and dispersion along the slit; filled circles are the data used by ROTCURVE while crosses were 
rejected by the program. The dashed curves show the intrinsic rotation and dispersion model profiles, and solid curves are the models after 
the seeing blur is applied. In panels (c) and (d) the solid curve overlies the dashed curve. The fitted rotation velocities and dispersions are 
Vrot = 149, 97, 16, 14 km s" 1 and 

CT 2d = 5, 5, 39, 30 km s 1 respectively. 



log and the 3" diameter magnitudes from the Capak 
catalog. For all but the largest, lowest-redshift galax- 
ies, aperture corrections are small, so we omit them. Of 
the 1440 TKRS galaxies, 1423 have matches with opti- 
cal magnitudes in the Capak et al. catalog. Fourteen of 
the remainder have magnitudes from the GOODS ACS 
catalog, and only 3 have neither. 

The Capak et al. HK' data are shallower than the op- 
tical data and the depth is not constant over the GOODS 
field. We used only aperture HK' magnitudes that have 
errors < 0.3 mag, yielding HK' for 919 TKRS galaxies. 
With this error cut, the fraction of TKRS galaxies with 
HK' magnitudes reaches 50% at HK' = 22.3 (AB). 

Our primary interests are to measure consistent rest- 



frame absolute blue magnitude Mb and color U—B, and 
a rest infrared magnitude Mj and R — J color. We con- 
vert to rest B and J because the observed filters cover 
or come relatively close to these restframe bands over 
< z < 1.5, minimizing extrapolation. We used the 
template-fitting AT-correction procedure described in pre- 
vious DEEP papers (Weiner et al. 2005; Willmer et al. 
2006) to convert observed colors into if -correction and 
restframe color. Briefly, this method uses rcdshiftcd tem- 
plates: 34 spectra of local galaxies with wavelength cov- 
erage 1200A - 1 [im (Kinney et al. 1996) for Mb and 
U - B; and 9 templates generated by PEGASE (Fioc 
& Rocca-Volmerange 1997) representing E, Sb and Sd 
SEDs at 1,5, and 10 Gyr for Mj and R— J. To analyze, 
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for example, a galaxy at z = 0.8 with R — I = 1.2, the 
templates are redshifted to z = 0.8 and their observed 
R— I colors synthesized. Since the templates' restframe 
U — B colors and Zf-corrections from I to B are known, 
we fit a low order polynomial to U — B as a function of 
R — I and evaluate this fit at R — I = 1.2 to find U — B 
for the galaxy under analysis. 

This procedure works best when the observed filter 
pair transforms relatively closely to the wavelength of 
the restframe filters, as is the case for R — I at z = 0.8 
into rest U — B. When the transformation is not close, 
there is larger scatter about the fit of U — B as a function 
of observed color, so choosing an appropriate filter pair 
is important. 

To derive the restframe U — B and Mb presented here, 
we used the filter pairs of Capak B — V,R — I, and I — z, 
switching over from one pair to the next at z = 0.6, 1.1 
respectively. For Az = 0.1 around the switchover red- 
shift we interpolate to make a smooth transition. For 
rest R — J and Mj we use the observed I — HK 1 filter 
pair and the HK' magnitude. 

In practice, the if-corrections are fairly stable. How- 
ever, certain filter pairs produce discordant results for 
restframe colors: e.g. ACS B — V and V — i produce col- 
ors for galaxies at z — 0.5 that are ~ 0.1 mag offset. This 
problem is most obvious with the ACS filter pairs and 
with the Capak V — R color. In general, a warning sign 
of the problem shows up as template mismatch in a color- 
color plot: for certain combinations of redshift and color 
pairs, galaxies and redshifted templates do not fall on ex- 
actly the same locus. The reasons for these mismatches 
are not clear; they may be a combination of small ze- 
ropoint shifts and imperfect filter/instrument response 
curves. The effect is small on absolute magnitudes, and 
relative color measurements at a given redshift are reli- 
able, but measurements of absolute color evolution over 
a large range in redshift can be affected. 

We found that using the Capak B — V, R — I, I — 
z color sets minimizes the effect of these mismatches; 
using the ACS filters produces slightly abnormal trends 
in restframe color with redshift. Comparing the absolute 
magnitudes and colors from the Capak color sets and 
from the fully independent ACS B — V, V — i, i — z 
sets, we find that 95% of the objects have a magnitude 
difference < 0.33 and a color difference < 0.26, yielding 
la errors of 0.12 in Mb and 0.09 in U — B. The errors 
on observed optical magnitudes from Capak et al. (2004) 
are 0.05-0.07 mag at the TKRS limit; additional scatter is 
propagated into the absolute magnitudes by the observed 
color error and the template scatter in the K-correction 
procedure. 

The median error on HK' magnitude for galaxies with 
kinematic measurements is 0.14 mag. Although we are 
using model rather than empirical SEDs to K-correct the 
IR magnitude, galaxy SEDs have less variety in the IR 
than the optical and so the K-correction depends only 
weakly on the specific templates used. The error on the 
IR observations dominates over the IR template scatter. 

3. SAMPLE MAGNITUDE, COLOR AND COMPLETENESS 

Of the 1440 TKRS galaxies, 1423 have magnitudes 
from Capak et al. (2004) and 14 have magnitudes from 
ACS (Giavalisco et al. 2004). 893 galaxies have magni- 
tudes and a linewidth determination, and an additional 



196 are kinematically unresolved. Of these 1089, 913 are 
brighter than Mb = —18, a limit we impose for Tully- 
Fisher fitting in Paper II. 681 of the 1089 with kinemat- 
ics have a reliable Mj and 647 of these are brighter than 
Mj = —19. Table [3 lists TKRS galaxies, positions, red- 
shifts, observed magnitudes, restframe magnitudes and 
colors, and linewidths. 

Figure 0] shows the restframe optical color-magnitude 
distribution of TKRS galaxies and the location of the 
1089 galaxies with and 348 without measured linewidths. 
The TKRS galaxies exhibit the color bimodality now 
well-established locally and to z ~ 1 (Strateva et al. 
2001; Bell et al. 2004; Wirth et al. 2004; Weiner et al. 
2005). Both red and blue galaxies show a shallow color- 
magnitude relation: brighter galaxies are redder in the 
mean. The majority of galaxies for which we fail to get 
a good linewidth are on the red side of the bimodality, 
U — B > 0.95, because those galaxies have weak or no 
emission lines (e.g. Weiner et al. 2005). Among the blue 
galaxies, linewidth failures are rare and fairly randomly 
distributed. We examined each failure and found that 
most are due to low line flux or night sky residuals. The 
failure rate in blue galaxies is small and there is no evi- 
dence that failures occur preferentially at large or small 
dispersions. 

At low redshift, blue galaxies in the TKRS sample show 
the well-known color-magnitude relation; fainter galax- 
ies are bluer (upper left panel of Figure 0J . At higher 
redshift, the faint limit is due to the TKRS magnitude 
selection; the selection limit in apparent R magnitude 
corresponds to a tilted line in the restframe Mb, U — B 
plane (e.g. lower left panel of Figure |3J. The tilt of this 
line changes with redshift. In Paper II we test the ef- 
fect of the selection limit on Tully-Fisher relations by 
defining an approximately mass-matched sample, apply- 
ing the z ~ 0.9 tilt at all redshifts and a magnitude cut 
that evolves similarly to I*. This cut is shown by the 
diagonal lines in Figure 0] 

The redshift completeness of the TKRS is ~ 70% to 
Rab = 24.4, discussed by Wirth et al. (2004). Redshift 
failures in the TKRS are mostly faint blue galaxies that 
are probably at z > 1.5. Here we discuss the linewidth 
completeness, i.e. the fraction of galaxies with good 
linewidths over galaxies with redshifts. 

Figure [5] shows the linewidth completeness as a func- 
tion of rest U — B color. The completeness is roughly 
constant for blue galaxies but drops sharply for red galax- 
ies. Because red galaxies are bright, a histogram of com- 
pleteness as a function of rest Mb for the whole sample 
actually drops at bright magnitudes. 

Red galaxies and blue galaxies are two distinct popu- 
lations and the bulk of our measurements are blue galax- 
ies, so in this paper we are constraining the behavior 
of the blue galaxy population. Figures El and show 
the linewidth completeness of the sample restricted to 
TKRS galaxies on the blue side of the color bimodality, 
U — B < 0.95. Linewidth completeness is nearly con- 
stant with either absolute or apparent magnitude, with 
a small tendency to be lower for galaxies bright in ab- 
solute magnitude (lower emission EW) or faint in ap- 
parent magnitude (lower apparent flux). The lack of 
any strong dependence on magnitude reflects the fact 
that intrinsically-faint galaxies tend to have high emis- 
sion equivalent width. 
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TABLE 2 

Catalog of TKRS galaxies and linewidth measurements 
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Note. — The complete version of this tabic is in the electronic edition of the Journal. The printed edition contains only a sample. Undefined 
values arc given as < —99. 



1 Redshift. 

^Apparent B RI (H K' ) AB aperture magnitudes (Capak ct al. 2004). Galaxies with only BI use GOODS ACS aperture magnitudes (Giavalisco 
et al. 2004). 

3 Restframe Mb, U — B, Mj and R — J AB magnitudes from the SED K correction procedure described in text. 

4 Log linewidth in km s _1 from integrated emission and its error. Galaxies assigned a log linewidth of 1.4 and error of 0.2 arc kincmatically 
unresolved, see text. 



Because completeness of the linewidth sample depends 
strongly only on color, we have a fair sample of the blue 
galaxy population, and our measurements will not be 
strongly affected by selection effects. Of course, any evo- 
lution measurement or interpretation will apply only to 
blue galaxies. 

4. PROPERTIES OF KINEMATIC MEASURES 

The limited spatial resolution of high-redshift spectra 
makes the interpretation of kinematic measures less ob- 
vious than in local galaxies. Our goal in this section is 
to establish the properties of our sample and determine 
the effect of observational limitations on the measures of 
dispersion and rotation. 

There has been some controversy over how well mea- 
sures such as integrated linewidth probe galaxy kinemat- 
ics. Modeling simulated observations of disk velocity 
fields with true circular velocity V c , Rix et al. (1997) 
found that for their observational parameters, the aver- 
age integrated linewidth < a > /V c = 0.6 with substan- 
tial scatter (<~ 0.15) depending on the unknown position 
angles and inclinations, where < a > incorporates an 
average over inclination but V c is the true value. For 
random orientations, < sin i >= 0.79, leaving an aver- 
age factor of 0.76 for <7/(VcSin i). From observations of 
local galaxies, Kobulnicky & Gebhardt (2000) found that 
the [O n] emission linewidth tracks the H I linewidth W20 
well for most galaxies, using a = O.28W20 (both uncor- 
rected for inclination), but the [O 11] linewidth is low in 
2/22 extreme cases. (The H 1 width W20 is the full- width 
of the H 1 profile at the 20% flux level.) This relation 
is roughly consistent with the Rix et al. relation since 



W 20 ^ 2(V c sin i + 15 km s" 1 ) (Tully & Fouque 1985). 

Lehnert & Heckman (1996) showed a tendency for 
edge-on IR-luminous starburst galaxies to have nuclear 
emission linewidth that is low compared to the expected 
velocity dispersion of the full potential; these galaxies 
are mostly fairly large and the nuclear emission does not 
probe the full potential. Barton & van Zee (2001) ob- 
served four blue compact dwarf galaxies and found that 
cr(Ha) falls below optical V rot and H 1 W50. However, 
their four galaxies average (T/V2.2 = 0.7, and the two 
with 21 cm data have cr/W^o = 0.27 (all quantities un- 
corrected for inclination). These arc actually fairly con- 
sistent with the < a > /V c offset of Rix et al. and the 
linewidth / H 1 width factor of Kobulnicky & Gebhardt. 
Rotation curves can also be affected: kinematic distor- 
tions and truncations could cause galaxies to have low 
measured rotation (Barton ct al. 2001; Kannappan & 
Barton 2004). 

Perhaps the most serious discrepancy between optical 
and H I width is found in a sample of 10 local blue com- 
pact galaxies (BCGs; Pisano et al. 2001), which show 
W2o(Hf3)/W 2 o(H 1) = 0.66±0.16. In terms of dispersion 
this implies <r(Hl3)/W2o(H 1) = 0.18. These galaxies are 
a fairly special population with high emission equivalent 
widths and small sizes, r e = 0.6 — 1.9 kpc, smaller than 
average for intermediate- z BCGs. They are most similar 
to NGC 4449, the furthest outlier in the sample of Kob- 
ulnicky & Gebhardt (2000), which has a much larger H 1 
extent than optical; and to the compact narrow emission- 
line galaxies (CNELGs, Koo et al. 1995) in their sizes. 
In the sample of Pisano et al. , the optical linewidths are 
most affected when I / l 7 2o(H 1) < 150 km s _1 , which corrc- 
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M„ AB, 0.8<z<l.l M B AB, l.Kz<1.61 



Fig. 4. — Absolute magnitude and color for TKRS galaxies in 4 redshift ranges. Blue Xes: galaxies with good linewidths; red circles: 
galaxies with poorly measured linewidth. The bulk of linewidth failures are in the red side of the color bimodality. The horizontal and 
diagonal lines show a magnitude and color selection that constructs samples whose depth is matched to the selection limit at z = 0.9 and 
an evolving L* . 



sponds to an unaffected optical a < 45 km s _1 ; smallest 
galaxies can be the most seriously affected. 

These types of underestimation of kinematic width 
may happen in high-redshift samples. They might ex- 
plain some of the discrepancies between previous TF evo- 
lution measurements, as Pisano et al. (2001) suggest, es- 
pecially in samples with high line EW. However, even at 
high redshift, the extreme BCGs or CNELGs are a mi- 
nority of blue galaxies (Koo et al. 1995). These galax- 
ies are quite small; even with 1-2 mag brightening, the 
Pisano et al. (2001) sample would be among the faintest 
of the z ~ 1 galaxies in the TKRS sample. The possibil- 
ity of BCG-induced offsets does argue in favor of using 
large samples which are not selected on emission line EW, 
and in which evolution can be measured internal to the 
sample. 

Ideal measures of high-redshift kinematics would yield 
2-d maps with high spatial resolution. This may soon 
be possible with adaptive optics and integral field spec- 
troscopy, but will remain impractical for large samples 



for some time. Here we carry out empirical tests on our 
sample and simulations of seeing-blurred observations to 
test the properties of linewidth measurements. 

4.1. Spatial extent of emission 

It is frequently suggested that high-redshift galaxies 
could have enhanced emission from their centers, such as 
a nuclear starburst, and that this could produce smaller 
linewidths because the emission does not probe the full 
potential, or comes from a part of the rotation curve 
that is rising (e.g. Lehnert & Heckman 1996). Simard 
& Pritchet (1998) found that ~ 25% of their galaxies 
had unresolved emission in their CFHT spectra. Rix et 
al. (1997) suggested that ~ 20% of their sample might 
have nuclear emission from high-density gas, based on 
the [O n] doublet ratio. However, both of these samples 
are selected to have high EW or blue colors, and are rela- 
tively low S /N compared to the TKRS spectra. Another 
possibility is that line emission could be dominated by 
one or a few large HII regions and so not probe the full 
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Fig. 5. — Linewidth completeness as a function of restframc 
U — B color. Solid line: galaxies with good lincwidths; dashed line: 
galaxies with poorly measured linewidth. The top panel shows 
absolute numbers and the bottom panel shows good or bad as a 
fraction of the total with rcdshifts. 
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apparent I mag AB, blue galaxies U-B<0.95 



Fig. 7. — Linewidth completeness for blue galaxies as a func- 
tion of apparent magnitude Iab- Solid line: galaxies with good 
linewidth; dashed line: galaxies with poor linewidth measure. The 
top panel shows absolute numbers and the bottom panel shows 
good or bad as a fraction of the total of blue galaxies with red- 
shifts. 
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Fig. 6. — Linewidth completeness for blue galaxies as a func- 
tion of absolute magnitude Mb- Solid line: galaxies with good 
lincwidths; dashed line: galaxies with poorly measured linewidth. 
The top panel shows absolute numbers and the bottom panel shows 
good or bad as a fraction of the total of blue galaxies with redshifts. 
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continuum FWHM, ", 0,8<z<l.l continuum FWHM, ", l.Kz<1.61 

Fig. 8. — Spatial extent along the slit for TKRS galaxies in 
four redshift ranges. The fitted gaussian FWHMs of emission and 
continuum are plotted. Filled circles have error on both FWHMs 
< 0.1", Small Xes have error on one or both FWHMs > 0.1", 
indicating an unreliable fit. The sizes are seeing-limited around 
FWHM ~ 0.7", and above that are consistent with equal sizes for 
emission and continuum. 



velocity field of the disk. 
We tested the spatial extent of emission in the TKRS 
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linewidth sample by fitting gaussian profiles along the 
slit. We collapsed the 2-d data in the wavelength di- 
rection over a 15 A range for emission and two 100 A 
ranges on either side of the line for continuum. We then 
fit gaussians using a non-linear least squares routine. 

The fitted spatial widths of emission versus continuum 
are shown in Figure 03 divided into four redshift ranges. 
The sizes are seeing-limited around FWHM ~ 0.7" and 
track each other well at larger sizes. The small mask- 
to-mask seeing variations in the TKRS data will tend 
to spread galaxies along the 1:1 line, but the relation 
continues far beyond FWHM ~ 1.0". Outliers from the 
1:1 line generally have a large error in one size measure- 
ment. While this comparison is somewhat limited by see- 
ing, there is no evidence for a large population of nuclear 
emission or single-HII-region emission. In fact there is 
strong evidence that emission and continuum size track 
each other well. 

4.2. What integrated linewidths measure: simulation of 
observations 

In order to understand what an integrated linewidth 
is really measuring, we simulated observations of a disk 
galaxy at z = 1, using an actual 2-dimensional intensity 
and velocity field from Fabry-Perot observations of the 
local galaxy NGC 7171. A similar project was carried 
out, also using Fabry-Perot velocity fields as input, by 
Rix et al. (1997). NGC 7171 has a fairly typical disk 
velocity field with mild spiral streaming motions, and 
abundant Ha emission, yielding a 2-d velocity field that 
samples the disk well with high resolution. It was ob- 
served in the Ha line with the Rutgers Fabry-Perot at 
the CTIO 4-m Blanco telescope, on October 27 and 28, 
1989 by T.B. Williams and R.A. Schommer 9 . The data 
were reduced into intensity, velocity, and dispersion fields 
by BJW using procedures described elsewhere (Palunas 
and Williams 2000; Weiner et al. 2001) and have been 
previously studied by Palunas (1996). NGC 7171 has a 
heliocentric velocity of 2740 km s~ x and we assumed a 
distance of 34 Mpc; its magnitude is Mb = —20.0. The 
resolution of these 2-d maps is 1.5", the inclination of 
NGC 7171 is 55°, and its rotation velocity corrected for 
inclination is V c = 189 km s _1 . It has a flat rotation 
curve and lies on the Tully-Fishcr relation of Sakai et al. 
(2000). 

The original Ha intensity, velocity and dispersion fields 
of NGC 7171 are shown in the left column of FigureEl the 
dispersion field is relatively featureless since the disper- 
sion is usually < 20 km s _1 except in the nucleus, where 
it is increased by seeing blur. We used these fields to re- 
construct a 3-d spectral cube. To simulate an observation 
at z = 1, we rescaled to the angular size at z = 1 , blurred 
each plane of the cube by a 0.7" gaussian to represent 
seeing, and resampled the cube onto DEIMOS pixels of 
0.118". We then refit through the cube to produce maps 
of restframe intensity, velocity, and dispersion. The right 
column of Figure El shows the resulting z = 1 maps that 
correspond to the unblurred maps of the left column. 

The blurred velocity and dispersion fields of Figure [5] 
illustrate an effect we alluded to in the discussion of ro- 
tation curve fitting in Section [2.3. 21 The beam smearing 

9 Cerro Tololo Inter- American Observatory is operated by 
AURA under contract to the NSF. 



due to seeing smooths out the rotation velocity gradient, 
but it produces a strong peak in the velocity dispersion 
in the center; where the unblurred velocity gradient is 
strong, the seeing mixes gas at different velocities to- 
gether. This peak in velocity dispersion is at high Ha 
intensity and so carries a high weight in integrated mea- 
surements. 

We used blurred velocity cubes to simulate observa- 
tions of integrated velocity dispersion over a range of 
galaxy inclinations and slit-galaxy relative position an- 
gle misalignments. To simulate a range of inclinations, 
we stretched the input NGC 7171 fields and rescaled the 
velocities before blurring and resampling. We did not al- 
low the dispersion to fall below 8 km s _1 , a typical sound 
speed for H II regions and gas in the ISM. We did not 
add noise or scale the flux to a real galaxy, as this is 
not a full simulation of real observations, but a test of 
the relative contribution of emission at different veloci- 
ties. Adding noise reduces the detectability of spatially 
resolved emission at low intensities but does not greatly 
affect the velocity distribution of the integrated emission. 

We then laid down typical 1" wide slits over a range 
of relative PAs and extracted a spectrum in a 2.3" (1.5 
FWHM) extraction window, using the intensities, veloc- 
ities and dispersions of pixels within the slit and window 
to build up an emission spectrum. The resulting velocity 
profile can be somewhat non-gaussian and flat-topped, 
especially when the slit is aligned with the galaxy major 
axis, but once the profile is convolved with the spec- 
trograph resolution, the non-gaussianity is small. We 
computed moments to determine the velocity dispersion 
corresponding to aid- Changing the extraction diameter 
within reasonable limits had little effect on the results, 
because the bulk of the signal comes from the central 
emission peak. 

FigureHHIplots the ratio of integrated linewidth to pro- 
jected circular velocity, a\d/(V c sin i), as a function of 
inclination and relative slit position angle. This ratio 
removes the sin i effect, isolating the properties of the 
integrated linewidth. There is a small residual effect of 
inclination, but the variation among the tracks is quite 
small. The exception is the nearly face-on i — 5° case, 
where the 8 km s _1 minimum dispersion contributes. 
The chances of observing a face-on, perfectly flat thin 
disk are very small, especially in the high-z universe 
where many galaxies are peculiar or irregular. Apart 
from this special case, the ratio varies by ±15% over the 
range of PA and inclination. 

At a given inclination, the linewidth for misaligned slits 
is lower, but the linewidth for a perfectly misaligned slit 
(APA — 90) is only 20% lower than the linewidth for a 
perfectly aligned slit (APA = 0) . The reason is that the 
observed linewidth comes largely from the central peak 
in intensity and dispersion, seen in the lower panel of 
Figure El The slit position angle does have a large effect 
on the observed rotation gradient, of course. Both these 
effects are seen empirically in our data in Section 14.41 
The ratio aid/(V c sin i) ~ 0.6 — 0.65 on average. When 
the distribution of inclinations is accounted for, the ratio 
is fairly similar to the distribution in Figure 7 of Rix et al. 
(1997): o~±d/V c ~ 0.55 in the mean, with a non-gaussian 
scatter due to the tail of nearly face-on model galaxies 
with low aid- 

The upshot of these simulations is that the seeing is 
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Fig. 9. — Emission intensity, line-of-sight velocity, and dispersion fields for NGC 7171 as observed locally (left column) and as it would be 
seen at z = 1 with 0.7" seeing and 0.118" DEIMOS pixels (right column). The scale bar is 30" on the left and 1" on the right. The velocity 
color scale is ±200 km s _1 and the dispersion color scale is 0-100 km s" 1 . The seeing blur smooths out the rotation velocity gradient, but 
causes a central spike of dispersion where emission from gas at different velocities is superimposed. 



a powerful integrator, smoothing over the disk and pro- 
ducing dispersion at the location of velocity gradients, so 
that specific observational parameters like slit PA, width, 
and extraction diameter do not have a strong effect on 
the integrated lincwidth. However, the slit PA does have 
a strong effect on the observed spatially-resolved rota- 
tion curve. Similar conclusions were reached by Erb 



et al. (2004) . Inclination produces significant scatter in 
lincwidth/velocity ratio merely through the sin i factor. 

A limitation is that we have modeled only a thin disk 
galaxy in orderly and planar rotation. Our comparison 
of observed rotation velocity and dispersion in Section 
15.11 below shows that many galaxies are not dominated 
by orderly rotation. However, even if the kinematics are 
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Fig. 10. — Ratio of integrated linewidth to projected circular 
velocity, <r 1( j/(V c sin i), as a function of inclination and slit position 
angle for simulated observations of NGC 7171 at z = 1. At a given 
inclination, the linewidth for a perfectly misaligned slit (APA = 
90) is only 20% lower than the linewidth for a perfectly aligned 
slit (APA = 0). The ratio takes out the inclination dependence, 
showing that the effect of inclination other than through sin i is 
small. The exception is the i = 5° case, where the sound speed of 
the gas, 8 km s _1 , contributes to the observed linewidth. 



tion is only practical with HST imaging, and galaxies 
may have stronger deviations from an ideal thin circu- 
lar disk geometry, so the relation between cllipticity and 
inclination is less direct. As groundwork for the Tully- 
Fisher study of Paper II, here we model the effect on 
disk Tully-Fisher relations when inclination and extinc- 
tion corrections are omitted. 
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not orderly, the effect of the seeing as an integrator over 
the galaxy's emission still operates. The net result is that 
the integrated linewidth is actually more robust than the 
rotation gradient against both unusual internal kinemat- 
ics and the observational effect of slit PA offsets. 

If the star formation and emission is more centrally 
concentrated than in NGC 7171, the relative contribu- 
tion of high velocity gas on the flat part of the rotation 
curve is reduced. However, because the central regions 
do have a velocity gradient, the integrated linewidth re- 
mains fairly high, and this is a second order effect. As a 
test, we artificially boosted the emission intensity inside 
1 kpc radius - the central clump of emission - by a factor 
of 4 and reran the simulations. The integrated linewidth 
decreased by just 4.5 ± 1%, over the range of inclina- 
tion and slit PA. Centrally concentrated emission must 
overwhelm the rest of the galaxy to lower the linewidth 
significantly. 

4.3. The effect of inclination corrections and the lack 

thereof 

The modeling of a blurred NGC 7171 showed that the 
strongest effect on the observed integrated linewidth is 
the inclination factor of sin i. Classic Tully-Fisher stud- 
ies correlate velocity and luminosity after correcting for 
inclination and extinction, and are generally restricted to 
relatively high-inclination galaxies with clear disk geom- 
etry. In this study we do not restrict the sample. In large 
samples of high-redshift galaxies, correcting for inclina- 



FlG. 11. — A simulated Tully-Fisher relation illustrating the 
effect of omitting inclination and extinction corrections. The black 
triangles are drawn from a Tully-Fisher relation with inverse TF 
slope -0.125 dex/mag and intrinsic scatter 0.0375 dex. The solid 
diagonal line is a linear fit to these points. The open circles are 
the same galaxies "observed" at randomly assigned inclinations 
i. The decorrected quantities are velocities V ro t sin i with 0.084 
dex observational error added and magnitudes Mg + Ag (i) with 
0.11 dex error. The dashed diagonal line is a linear fit to the 
decorrected points. The lines at right indicate the decorrection 
tracks for a galaxy with log V ro t = 2.2 at i = 10°, 30°, 60°, 90° . 
Edge-on galaxies move above the baseline inverse-TF relation due 
to extinction; face-on galaxies fall below the baseline due to low 
projected velocities. For B band extinction, galaxies at the median 
inclination i = 60 are nearly unchanged. However, the scatter 
around the TF relation becomes large and non-gaussian. 



Figure ^2 shows a model Tully-Fisher relation with 
and without inclination and extinction corrections. The 
filled triangles are a simulated "true" Tully-Fisher re- 
lation with a realistic distribution of magnitudes (drawn 
from our sample in < z < 0.5), an inverse-TF log V(M) 
slope of —0.125 dex/mag, and an intrinsic scatter of 
0.0375 dex (0.3 mag). Each galaxy is assigned an in- 
clination drawn from a uniform distribution in cos i; 
we do not modify the inclination distribution to reflect 
the small bias against edge-on galaxies in a magnitude- 
selected sample, and surface brightness selection against 
face-on galaxies. The open circles are the decorrected 
values log V ro t,obs and MB,obs, with errors added from 
gaussian distributions of 0.084 dex and 0.11 mag. The 
extinction correction (to face-on) is Ab = 7s log(a/6), 
with 7s = 1.57 (Sakai et al. 2000, neglecting the lu- 
minosity dependence of extinction), and we assume the 
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galaxies have intrinsic thickness:radius 1:5. The veloc- 
ity is Vrot sin i with a random-motion component of 25 
km s _1 added in quadrature, which prevents the face-on 
galaxies' velocity from falling to near zero. The solid and 
dashed lines show the best linear fits respectively to the 
original and decorrected points with Mb < —18 (using 
a fit method that compensates for scatter, described in 
Paper II). The fitted slopes are -0.126 ± 0.003 for the 
original sample and —0.105 ± 0.021 for the decorrected 
sample. 

The lines at right show how decorrection moves 
a galaxy with log V ro t = 2.2 for inclinations 
10°, 30°, 60°, 90°. The undoing of these inclination and 
extinction corrections moves the galaxies onto tracks 
nearly parallel to the original TF relation e.g. all the 
i = 90° galaxies are shifted to fainter magnitudes, so 
appear above the original TF. For jb — 1-57 and a TF 
slope —0.125, the corrections cancel at i = 60.5°, which 
fortuitously is almost exactly the median inclination of 
i = 60° for a randomly oriented sample. Thus in i?-band, 
the effect of undoing corrections on the TF intercept is 
relatively small. For the galaxies on a single inclination 
track, the slope is unchanged, but for the entire sam- 
ple, the slope of the uncorrected inverse TF relation is 
shallower. 

As shown in Figure ^2 the scatter induced by lack of 
corrections is non-gaussian, although the non-gaussianity 
is blurred somewhat by the observational errors. The ef- 
fect on Tully-Fisher fitting depends on the magnitude 
distribution and limit of the sample. The fitted TF rela- 
tion changes due to the correlated anti-corrections scat- 
tering across the TF ridgeline, the non-gaussian scat- 
ter, and the elimination of cdge-ons by the magnitude 
limit. For the realization shown here, a fit to the galaxies 
brighter than Ms.obs = — 18 yields a velocity intercept 
shifted down by 0.11 dex at Ms.obs — —20, a slope of 
—0.105 ± 0.021, and a greatly increased scatter of 0.19 
dex. We emphasize that these values only apply for the 
idealized thin disk assumptions, which likely do not apply 
to many galaxies in the TKRS - for example, elongated 
but non- rotating galaxies like TKRS 5627 and 10023, 
shown in Figures ^ an d [21 

Local Tully-Fisher samples require a luminosity or 
velocity-dependent extinction correction; extinction is 
higher in larger galaxies (Giovanelli et al. 1995; Tully 
et al. 1998). Undoing this correction means that an un- 
corrected inverse TF sample has a V{M) slope steeper 
than the true slope. In the B band the effect is about 10- 
15% of the magnitude range of the sample (Tully et al. 
1998), so the uncorrected inverse TF slope is ~ 10— 15% 
higher. This could counteract the slope-shallowing ef- 
fect described above. However, one must tread carefully 
when applying local extinction corrections to distant 
samples. Both the extinction effect and the inclination- 
induced offsets are reasons to attempt comparisons in- 
ternal to the sample, spanning a range of redshifts. 

The general effects of the lack of correction are small 
shifts in intercept and slope, and a substantial increase 
in scatter. If the TKRS galaxies with linewidths were all 
close to the thin rotating disk model, applying inclina- 
tion and extinction corrections would greatly reduce the 
scatter in the Tully-Fisher relations fitted in Paper II. In 
fact, in Paper II we find that applying either the sin i 
inclination correction, or both inclination and extinction 



corrections, does not reduce the scatter in Tully-Fisher 
relations. A small amount of this can be due to statisti- 
cal errors on inclination, but most is due to the diversity 
of kinematic properties, galaxies with non-disk shapes, 
and motions that are not strictly rotation, as we show in 
Section El 

4.4. Empirical dependence on slit position angle 

Rotation curves of inclined disks should be harder to 
measure when the slit is not aligned close to the major 
axis. The modeling of NGC 7171 suggests that integrated 
velocity linewidths are not as severely affected by slit 
misalignments. The effect can be tested empirically for 
the TKRS sample since HST imaging allows measures of 
the galaxy position angles. 

We used a compilation of measurements from ellipse 
fitting to the ACS i-band galaxy images, yielding ellip- 
ticity, position angle, and half-light radius. Segmenta- 
tion images output by the SExtractor program (Bertin & 
Arnouts 1996) were used to mask out neighboring galax- 
ies and the IRAF task ellipse was run to fit elliptical 
isophotes. 10 Further results on the galaxy sizes will be 
reported separately (Melbourne et al. 2006). Here we 
use the ellipticities and position angles to determine slit 
misalignments. 
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Fig. 12. — Integrated linewidth log <r 1(J as a function of position 
angle misalignment: offset between slit PA and major axis PA 
from ellipse fits to HST/ACS images. Large points show the mean 
and RMS in bins of angle. The small points at log <T lt j = 1.4 are 
the kinematically unresolved objects. There is no visible trend of 
measured 1-d linewidth with PA offset. 



FigurelT2lshows integrated linewidth log <r 1( i as a func- 
tion of the misalignment between slit and image position 
angle. The large points show the mean and RMS in bins 

10 IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under a cooperative agreement with the National Science 
Foundation. 
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of PA offset. (The slitmasks were designed based on PAs 
from ground based imaging before the HST catalog was 
available, so the slits are not exactly aligned with HST 
PAs.) Slit misalignment has no visible effect on linewidth 
measured from the 1-d spectrum. A fraction of the galax- 
ies in this plot are nearly round, so that the PA is not 
very meaningful; removing these objects does not change 
the conclusion. A small trend with PA offset might be 
expected from the models in Figure^B but galaxies that 
are not orderly rotating disks, and galaxies whose ellip- 
ticity and kinematics are misaligned, will weaken such a 
trend. 



TABLE 3 

Catalog of TKRS rotation and dispersion fits 



2 r 



40 60 
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Fig. 13. — Rotation velocity log V ro t as a function of position 
angle misalignment, offset between slit PA and HST/ACS major 
axis PA. Round galaxies with ellipticity e < 0.25 are plotted as 
Xes, non-round galaxies are circles. Large points show the mean 
and RMS of non-round galaxies in bins of angle. Extracted rot ation 
velocities are higher for aligned slits; there is a deficit of high V ro t 
for misaligned slits. 



For the subsample of galaxies with 2-d rotation fits, 
Figure IT31 shows the rotation velocity log V ro t as a func- 
tion of the slit misalignment angle. Again, large points 
show the mean and RMS in bins of PA offset. Here, there 
is a strong effect: measured rotation velocity declines as 
misalignment increases, and there is a lack of high V ro t 
objects for misaligned slits, as expected. However, note 
that even aligned slits have a fair number of low V ro t 
objects. 

Clearly, measuring rotation velocity from a mis- 
aligned slit is subject to large error, especially since 
ROTCURVE does not use a full 2-d model. For pur- 
poses of studying the properties of rotation velocity 
and dispersion in the next section, we split the sample 
into "aligned" and "misaligned" with |APA| < 40° and 
> 40°. We also exclude the nearly-round galaxies with 
ellipticity e < 0.25 from the "aligned" subsample. 

5. PROPERTIES OF SPATIALLY RESOLVED KINEMATICS 
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0.6 
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0.4855 


2.1 


110.2 
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edge-on 
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disk 



Note. — The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample. 

1 Redshift. 

^Absolute value of the offset between DEIMOS slit position angle and 
photometric major axis PA in HST/ACS / images, degrees. 

^ Line-of-sight rotation velocity from seeing-compensated modeling of 
the 2-d spectra, kms -1 . For |APA| > 40°, rotation velocities are 
strongly affected by slit misalignment. 

^ Line-of-sight dispersion from seeing-compensated modeling of the 2- 
d spectra, km s~ 1 . Values are near-quantized in units of 5 km s _1 due 
to the model grid. 



Morphological type from visual classification. Only galaxies with 
|APA| < 40 and ellipticity > 0.25 were classified. 



5.1. Rotation versus dispersion dominated galaxies 

The two galaxies whose rotation and dispersion profiles 
are shown in Figure [3] exemplify a trend we found in the 
entire sample. In order to fit the kinematics, we had to 
include both velocity V ro t and dispersion cr 2t i in the 2- 
d modeling of ROTCURVE, discussed in Section 
Quite a few galaxies like TKRS 5627 and 10023 show a 
low rotation velocity but a significant dispersion. 

Table |3 lists galaxy ID, redshift, slit PA offset, and the 
fitted line-of-sight rotation and dispersion for the 380 
galaxies with good ROTCURVE fits. It also includes 
morphological classifications for a subset of the galaxies, 
discussed further in Sectional Because there are galaxies 
which appear not to be consistent with inclined thin ro- 
tating disks, we report the line-of-sight velocities rather 
than applying an ellipticity-based inclination correction. 
If a significant number of high-redshift galaxies are not 
circular rotating disks, the chain of reasoning from axis 
ratio to inclination to sin i correction breaks down. 

The simulated seeing-blurred velocity fields in Section 
14. 21 suggest that observing a rotating galaxy with the slit 
along the minor axis yields a low velocity gradient but a 
central peak in velocity dispersion. These are found in 
our sample, but TKRS 5627 and 10023, shown in Figures 
and|3 are different. They exemplify a class of galax- 
ies which have low rotation, high dispersion, and aligned 
slits (|APA| < 40°), so that the lack of rotation is not 
due to slit misalignment. These galaxies also frequently 
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show dispersion profiles which are not sharply peaked, 
although the limited spectral resolution of TKRS makes 
dispersion profiles rather noisy. 

We classified galaxies into rotation- and dispersion- 
dominated, based on whether V ro t > aid or vice versa. 
Figure O plots aid versus V rot ; the placement of the 
dividing line is somewhat arbitrary, but intended to sep- 
arate galaxies whose main kinematic support is from or- 
derly rotation, and galaxies where dispersion (or the dis- 
ordered gas motions that aid may represent) plays a ma- 
jor role. When roundish galaxies with e < 0.25 and those 
with misaligned slits, |AP^4| > 40°, are excluded, there 
are 185 galaxies with good 2-d fits from ROTCURVE. 
Of these, 118 are rotation dominated and 67 (36%) are 
dispersion dominated. Thus two-thirds of the subsample 
have kinematics basically similar to local rotating galax- 
ies - which was by no means guaranteed - but over one- 
third of this subsample have kinematics that are incon- 
sistent with a thin rotating disk model. 

The nature of the dispersion dominated galaxies is a 
puzzle. They are not literally dispersion-supported ellip- 
tical galaxies. Comparing rotation and dispersion dom- 
inated galaxies (RDGs and DDGs) in magnitude and 
color, the DDGs are only slightly different, 0.05 bluer 
in median U — B and equal in median brightness, and 
their morphologies are late type; we discuss DDG and 
RDG properties further in Section The DDGs' ve- 
locity dispersion is unlikely to indicate a literal dynami- 
cally hot pressure-supported tracer population like stars 
in spheroidals. We are measuring velocities from nebular 
lines that come from ~ 10 4 K gas, so the dispersion in 
individual H n regions is only 8-10 km s . The larger 
values of aid found in DDGs could come from the relative 
motions of the ensemble of gas clouds - chaotic, disor- 
dered, or out-of-plane velocity fields that are smoothed 
over by the seeing, or that have strong misalignments 
between photometry and kinematics. However, the ge- 
ometry of the motions and whether they will eventually 
dissipate and settle into a disk cannot be resolved in these 
spectra. 

If the seeing blur caused a bias against the detection 
of rotation in small objects, it could make objects ap- 
pear dispersion dominated by masking the true rotation 
gradient. In principle, the ROTCURVE fitting method 
is seeing-compensated. For a very small object with a 
rotation gradient that suffers greatly from seeing blur, 
it yields a shallow minimum in x 2 ; models with an in- 
trinsic rotation gradient and with a dispersion both fit 
the data, so there is a large error on V ro t rather than a 
bias. In practice, we rejected objects with a very small 
diameter of velocity data, < 1". In Sectional we discuss 
the sizes of RDGs and DDGs and show that the DDGs 
are only mildly smaller, so they are not simply caused 
by an observational bias against rotation. Intrinsic non- 
roundness or errors in the ellipticity and PA measure- 
ments could cause a few face-on disks to appear in the 
DDG sample, but there are too many elongated DDGs 
to explain away. 

Dispersion dominated galaxies occur over the full red- 
shift range of our sample and are not strictly a high- 
rcdshift population. Figure [TBI plots magnitude against 
redshift for the sample with 2-d ROTCURVE model- 
ing, with points coded for rotation or dispersion domina- 
tion. DDGs exist at low redshift, but are relatively faint 




log V rol , from rotation curve 



Fig. 14. — 2-d dispersion versus rotation velocity for galaxies with 
2-d ROTCURVE modeling. Open circles: dispersion dominated 
galaxies, V ro t < "m, wn -h e > 0.25. Filled triangles: rotation 
dominated galaxies, Vrot > o"2ti> with e > 0.25. Small Xes: round 
galaxies or misaligned slits. We divide the galaxies into rotation 
and dispersion dominated at V ro t = <J2d- The values of o^d are 
quantized in steps of 5 km s by the fitting method. 



there; it is possible that the dispersion/rotation fraction 
is higher at high redshift, although the sample size and 
selection effects make any conclusion preliminary. 

5.2. Comparison of integrated linewidths and spatially 
resolved kinematics 

A key question for using linewidths of integrated emis- 
sion as a characteristic velocity scale is how well they 
correlate with spatially resolved kinematic measures. 
In Section 14.21 we tested models of seeing-blurred in- 
tegrated emission from a rotating disk and found that 
linewidth aid correlated well with projected circular ve- 
locity V c sin i; the chief scatter between aid and true V c 
is induced by inclination. In this section we test the rela- 
tion of linewidth to rotation and dispersion from spatially 
resolved fits empirically using the 1-d and 2-d kinematic 
fits in TKRS. These quantities are line-of-sight, so no 
inclination correction is applied. 

The existence of rotation and dispersion dominated 
galaxies suggests that neither 2-d measure, Vrot nor aid, 
by itself adequately represents the velocity scale of every 
distant galaxy. Figures ITTjl and IT7I plot the 2-d velocity 
and dispersion, V ro t and aid, against the 1-d integrated 
linewidth, aid- Point types are coded by whether the 
galaxy's 2-d kinematics are rotation or dispersion dom- 
inated. The clustering of points with log aid = 1-4 are 
galaxies kinematically unresolved in 1-d. 

Figures [TBI and IT7I are the flip-side of each other. Each 
2-d quantity, when it dominates, correlates with the 1-d 
measure aid, but falls well below aid when it does not 
dominate. Either used alone would leave some galaxies 
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Fig. 15. — Magnitude versus redshift for galaxies with 2-d 
ROTCURVE modeling. Open circles: dispersion dominated galax- 
ies, Vrot < f2d) with e > 0.25. Filled triangles: rotation dominated 
galaxies, Vrot > om, with e > 0.25. Small Xes: round galaxies or 
misaligned slits. Dispersion dominated galaxies exist at low red- 
shift, but are relatively faint. There is not strong evidence for 
redshift evolution in the rotation/dispersion fraction, but possibly 
bright dispersion dominated galaxies are only found at high red- 
shift. 



Fig. 16.— Rotation velocity from the 2-d ROTCURVE mod- 
eling, log V r ot, versus integrated linewidth log cr ld . Open cir- 
cles: dispersion dominated galaxies, V ro t < &2d, with e > 0.25. 
Filled triangles: rotation dominated galaxies, Vrot > o" 2£ j, with 
e > 0.25. Small Xes: round galaxies or misaligned slits. The 
galaxies at log <r ld = 1.4 are kinematically unresolved in the 1- 
d linewidth. The diagonal lines are the 1:1 line and the Rix et 
al. (1997) a = 0.6V C line. For dispersion dominated galaxies, 
captures the velocity scale but V ro t goes to small values. 



with egregiously small velocities. 

For rotation dominated galaxies, there is a fair amount 
of scatter between V ro t and aid, but in general these 
quantities are correlated, with aid offset below V ro t- But 
for dispersion dominated galaxies, V ro t is low and does 
not reflect the kinematic support seen in aid- Con- 
versely, aid correlates well with aid in dispersion domi- 
nated galaxies. But rotation dominated galaxies can be 
fit without requiring 2-d dispersion, so aid falls below 
aid- In the RDGs, aid is capturing the seeing- induced 
dispersion. 

The inadequacy of each 2-d quantity V ro t and aid by 
itself motivated the design of ROTCURVE, which fits 
both simultaneously. To represent a characteristic veloc- 
ity derived from the 2-d ROTCURVE fits, we define a 
"combined velocity scale" Sk- 

Si = KV rot 2 + aid 2 , (5) 

where K is a constant < 1. The scale Sk combines rota- 
tion and pressure or random/disordered-motion support, 
just as Tully-Fisher studies have done to model an ob- 
served integrated H I linewidth (e.g. Tully & Fouquc 
1985). The simplest combination is Si.o, with K = 1. 
This is in some sense an estimator of rotation velocity 
plus a turbulent-motion correction. 

However, setting K < 1 makes Sk a better estima- 
tor of the velocity dispersion a of the entire potential, 
because rotation velocity and dispersion contribute dif- 
ferently to the support of a gravitationally bound sys- 



tem. For example, in an isothermal potential with a flat 
rotation curve, a population with isotropic velocity dis- 
persion and no net rotation has a = V c /V2, suggesting 
that with K = 0.5, So. 5 is a good estimator of a of the 
potential. In general, for a spherically symmetric tracer 
distribution with density oc r~" and isotropic velocity 
dispersion, if V c varies slowly with radius, a = V c /yfa 
(Binney & Tremaine 1987). Real galaxies in the range 
at or beyond the peak of the baryonic contribution to 
the rotation curve, where rotation curves are measured, 
have a = 2 — 3 for a spherized mass distribution, so 
K = 0.3 — 0.5 is reasonable. 

Figures ED and ED plot the combined velocity scales 
5*1.0 and So. 5 against 1-d linewidth aid- Combined ve- 
locity alleviates the shortcomings of V ro t and aid, elimi- 
nating artificially low velocity measurements. Both Si.o 
and So. 5 show strong correlations with aid- For the sub- 
sample with aligned slits and e > 0.25, the mean offset 
log Si.o — log fid is 0.15 dex with RMS 0.17 dex; for So. 5 
the mean and RMS log offset are 0.06 and 0.14. The ro- 
tation estimator Si.o shows its rotation velocity heritage; 
rotation dominated galaxies have Si.o > ^ld, relatively 
close to the a = 0.6V C line of Rix et al. (1997). RDGs 
and DDGs are offset in Si.o; the offset in So. 5 is smaller. 
When the dispersion-analog So. 5 is plotted against aid, 
the correlation with aid tightens, and the offset between 
So. 5 and aid is reduced, because both are estimating a 
velocity dispersion of the potential. 

The good correlation between 1-d integrated linewidth 
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log (J ld , km/s, 1-d dispersion 

Fig. 17.— Dispersion from the 2-d ROTCURVE modeling, log 
<*ldi versus integrated linewidth log <r l£ j. The G2d values are quan- 
tized by fitting on a 5 km s — 1 model grid. Open circles: dispersion 
dominated galaxies, Vrot < o"2di with e > 0.25. Filled triangles: 
rotation dominated galaxies, V ro t > &2d, with e > 0.25. Small Xes: 
round galaxies or misaligned slits. The diagonal line is the 1:1 line. 
For dispersion dominated galaxies, aid and ct2d are correlated; for 
rotation dominated galaxies, <t 1( j captures the velocity scale while 
&2d becomes small. 



and the combined velocity scale from 2-d rotation and 
dispersion fits has two major implications. First, 1-d 
linewidth is a reliable representation of the more sophisti- 
cated 2-d kinematics. As suggested by the simulations in 
Section Ol the seeing is a powerful integrator. The 1-d 
linewidth can be used as a kinematic measure when esti- 
mating dynamical masses, although the aid ~ 0.614 and 
the substantial scatter should be kept in mind. Second, 
2-d rotation velocity alone does not completely account 
for the kinematic support for all galaxies. For galaxies 
that are selected to have visibly rotation-dominated kine- 
matics, V ro t is adequate, but for broader samples, V ro t 
alone will not suffice. Constructing a Tully-Fisher rela- 
tion using V ro t only will push galaxies with U2d > V ro t 
to erroneously low velocity, as seen in Figure Hrjl 

6. DISCUSSION: THE NATURE OF ROTATION AND 
DISPERSION DOMINATED GALAXIES 

The existence of rotation and dispersion dominated 
emission line galaxies - essentially two modes of dynam- 
ical support, or a sort of kinematic morphology - raises 
the question of their nature and evolution. Rotation is 
well understood in local galaxies; dispersion dominated 
galaxies are the real puzzle. The size of the TKRS sample 
yields enough RDGs and DDGs to compare the physical 
properties and morphologies of these galaxies. 

Several other studies of intermediate and high redshift 
galaxies have found objects inconsistent with rotation, 
similar to the DDGs. Simard & Pritchet (1998) found 
25% of their galaxies had "kinematically anomalous" [O 




log (7 ld , km/s, 1-d dispersion 

Fig. 18. — Combined velocity, Si.o 2 = Vrot 2 + cr 2 d 2 : log Si.o 
versus integrated linewidth log o~id- Open circles: dispersion dom- 
inated galaxies, V ro t < om, with e > 0.25. Filled triangles: ro- 
tation dominated galaxies, V ro t > aid-, with e > 0.25. Small Xes: 
round galaxies and those with large error on V ro t and a 2d- The 
diagonal lines are the 1:1 line and the Rix et al. (1997) a = 0.6V C 
line. Linewidth and Si.o are correlated, falling closer to the 1:1 
line for dispersion dominated and to the Rix et al. line for rotation 
dominated galaxies. 



II] emission. These objects had centrally concentrated 
emission; we find below that DDGs are smaller in both 
continuum and emission. Erb et al. (2004) found that 
spectra of 13 elongated galaxies at z ~ 2 rarely show 
rotation. Bershady et al. (2005) obtained HST/STIS 
spectra of six luminous compact blue galaxies and found 
low V/a; the STIS velocity dispersion confirmed ground- 
based integrated linewidth. From integral field spec- 
troscopy, Flores et al. (2006) found that ~ 2/3 of 35 
galaxies at 0.4 < z < 0.75 had either mildly decentered 
or non-rotating kinematics; some of these may be dis- 
persion dominated. Our TKRS survey establishes that 
the dispersion dominated galaxies are relatively common 
and not restricted to specially selected samples. 

It is unclear what the local counterparts of disper- 
sion dominated galaxies could be. Many local studies 
of galaxy kinematics, in order to study the distance scale 
or rotation curves and dark matter, have selected against 
galaxies with disordered kinematics. However, even very 
small galaxies down to Mb ~ —15—16 are generally 
dominated by rotation (e.g. Swaters et al. 2002). Tully- 
Fisher studies have considered the contribution of tur- 
bulent motions to the linewidth; Tully & Fouque (1985) 
claim that turbulent motions are only significant for very 
small rotation velocities V < 25 km s _1 , but this is based 
on an idealized model. Kannappan & Barton (2004) find 
a larger percentage of kinematically anomalous galax- 
ies in close pairs; these mostly show distorted or trun- 
cated rotation curves, but a few are unusual enough that 
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Fig. 19. — Combined velocity, So.5 2 = 0.5V ro t 2 + cr 2d 2 : log 
So. 5 versus integrated linewidth log cri^. Open circles: dispersion 
dominated galaxies, V ro t < &2di with e > 0.25. Filled triangles: 
rotation dominated galaxies, V ro t > &2d, with e > 0.25. Small Xes: 
round galaxies and those with large error on V ro t and 02d- The 
diagonal lines are the 1:1 line and the Rix et al. (1997) a = 0.6V C 
line. Linewidth and So. 5 are correlated; the 0.5 prefactor makes the 
combined velocity width a better estimate of velocity dispersion, 
so that the correlation is tighter and the galaxies closer to the 1:1 
line. 



seeing-blur could make rotation gradients undetectable. 
It is also possible that evolution makes bright dispersion 
dominated emission line galaxies rare at the present day. 

6.1. Comparative Properties of Rotation and Dispersion 
Dominated Galaxies 

Rotation and dispersion dominated galaxies do exist 
over the full rcdshift range of our sample, as seen in Fig- 
ure ED There is some sign that at z < 0.4, DDGs are 
only found at faint magnitudes. This might reflect an 
evolution either in kinematics or in magnitude and emis- 
sion strength. To compare the properties of luminous 
RDGs and DDGs, and to avoid being pulled by faint ob- 
jects visible only at low redshift, we show distributions 
of galaxy properties for a sample restricted to z > 0.4, 
and with aligned slits and ellipticity e > 0.25 as before. 

The distributions of z > 0.4 RDG and DDG restframe 
magnitude and color Mb and U — B are shown in Fig- 
ures |201 and El At z > 0.4, the two types have in- 
distinguishable magnitude distributions. The DDGs are 
bluer in the median than the RDGs, but only by 0.05 
mag. Comparing the RDG and DDG distributions with 
the Kolmogorov-Smirnov statistic gives a 93% probabil- 
ity that they are identical in Mb and a 15% probability 
that they are identical in U — B. The DDG color dis- 
tribution covers almost the whole of the range of the 
blue galaxy population, so it appears that DDGs are not 
solely unusual blue objects, e.g. compact H II regions or 
extreme starbursts. 



Fig. 20. — Magnitude distribution in restframe Mg of rotation 
and dispersion dominated galaxies at z > 0.4, with aligned slits 
and e > 0.25. Solid line: rotation dominated galaxies, V ro t > &2d- 
Dashed line: dispersion dominated galaxies, Vrot < (J 2d- At high 
redshift, RDGs and DDGs are equally bright. 



rot dominated 
disp dominated 



U-B, z>0.4 



Fig. 21. — Color distribution in restframe U — B of rotation 
and dispersion dominated galaxies at z > 0.4, with aligned slits 
and e > 0.25. Solid line: rotation dominated galaxies, Vrot > cr 2 d- 
Dashed line: dispersion dominated galaxies, Vrot < &2d- RDGs 
are slightly redder than DDGs. 



The size distributions of z > 0.4 RDGs and DDGs are 
shown in continuum and in emission, in Figures l2"2l and 
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Fig. 22. — Radius distribution in half-light radius R^i, kpc, from 
HST images, for rotation and dispersion dominated galaxies at z > 
0.4, with aligned slits and e > 0.25. Solid line: rotation dominated 
galaxies, Vrot > &2d- Dashed line: dispersion dominated galaxies, 
Vrot < &2d- RDGs are larger in the mean than DDGs, though the 
size distributions overlap. 



Fig. 23. — Emission extent: the total range (diameter), in arcsec, 
of data used in the rotation/dispersion fit in the Keck spectra, 
for rotation and dispersion dominated galaxies at z > 0.4, with 
aligned slits and e > 0.25. Solid line: rotation dominated galaxies, 
V ro t > &2d- Dashed line: dispersion dominated galaxies, V ro t < 
d2d- RDGs are larger in the mean than DDGs, though less so than 
in half-light radius. The ranges in emission extent are similar, 
implying that DDGs are real and not just objects too small to 
detect rotation in ground-based spectra. 



1231 Figure plots the physical half-light radius Rhi in 
kpc from the HST ACS / images, while Figure l2~3l plots 
the observed emission diameter in arcsec, D em , using 
the range of data in the Keck spectra deemed good by 
ROTCURVE fitting. 

Here there is a difference between DDGs and RDGs. 
Dispersion dominated galaxies are smaller in both size 
measures. The offset is more noticeable in the physical 
radius Rhi- Although the range of sizes for DDGs is 
about as broad as for RDGs, the DDGs are somewhat 
smaller in the average: the medians of Rhi are 3.7 and 
4.7 kpc for DDGs and RDGs. The DDGs are similar in 
luminosity to RDGs but smaller in radius, so they are 
somewhat more compact. 

The emission diameter D em plotted in 1231 is the an- 
gular extent from the first to last row (inclusive) in the 
Keck spectrum of the data that ROTCURVE deemed 
good for fitting. There is a mild offset, with rotation 
dominated galaxies being larger in the average, but the 
RDGs and DDGs span the same range of angular sizes. 
The medians of D em are 2.26" and 2.38" for DDGs and 
RDGs. The distribution of D em is limited at the small 
end because the rotation curve sample was selected by 
enforcing a lower limit in ground-based angular size of 
FWHM = 0.94". The K-S statistic gives a 2% probabil- 
ity that the RDG and DDG distributions are identical 
in Rhi and a 12% probability that they are identical in 
D em . From the K-S tests, the only firm difference be- 
tween RDG and DDG properties is in the physical radius 
Rhi- The two types of galaxy are not as different in D em 
as they are in HST- measured physical radius; this could 
be caused by the selection on ground-based size; by see- 



ing blur; by differences in the falloff of the light profiles; 
and/or by threshold effects, since D em measures where 
the emission strength drops below some value that yields 
a reliable velocity. 

The chief point of Figure [231 is that the DDGs and 
RDGs in our rotation curve sample span a similar range 
of angular extent as measured on the 2-d spectra. This 
implies that the dispersion dominated galaxies are gen- 
uinely non-rotating in their mean emission. It argues 
against the possibility that DDGs are simply an obser- 
vational bias in small galaxies, caused by emission which 
does not extend far enough to detect rotation, alluded to 
in Section 15. II However, such a bias might re-appear in 
samples without an angular size lower limit. 

6.2. Morphologies of Rotation and Dispersion 
Dominated Galaxies 

To explore further any systematic differences between 
rotation and dispersion dominated galaxies, we made a 
morphological classification of the sample of 116 RDGs 
and 67 DDGs that have aligned slits and ellipticity 
> 0.25. We inspected ACS /-band images from the 
GOODS-N data. The classifications are cataloged in Ta- 
ble 

Because the physical resolution is limited and high red- 
shift galaxies are known to exhibit a greater degree of 
morphological peculiarity (e.g. Abraham et al. 1996), 
we defined a simplified system of morphological types. 
We found no spheroidals in the RDG/DDG sample. The 
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TABLE 4 

Morphologies of rotation and dispersion dominated 
galaxies 



Type 


RDG 


RDG 


DDG 


DDG 




all z 


z > 0.4 


all z 


z > 0.4 


spheroidal 


0% 


0% 


0% 


0% 


disk 


42% 


45% 


22% 


27% 


edge-on disk 


16% 


8% 


2% 


0% 


irregular 


31% 


32% 


43% 


41% 


chain/hyphen 


7% 


10% 


25% 


22% 


merger 


3% 


5% 


8% 


10% 


total 


116 


74 


67 


49 



classes are: disks, Sb-Sd; irregulars, both low and high 
surface brightness; edge-on disks; "chain" and "hyphen" 
galaxies; and mergers. Chain galaxies are elongated ob- 
jects with bright knots (e.g. Cowie, Hu, & Songaila 
1995). Hyphen galaxies are thin, small but elongated, 
with little substructure. We only defined galaxies as 
mergers if they clearly had multiple components, tidal 
tails, or major disturbances associated with mergers. 
Many high-z galaxies are somewhat peculiar or asym- 
metric and we did not count these as evidence of merg- 
ers. The percentages of morphological types for RDGs 
and DDGs are listed in Table 0J for both the full sample 
and the sample restricted to z > 0.4, where RDGs and 
DDGs have similar luminosity. 

The clearest results from Table 0] are that rotation 
dominated galaxies are generally disky, or irregular, and 
very rarely chain or hyphen galaxies. Dispersion dom- 
inated galaxies sometimes appear disky, but are more 
often irregular, and a substantial number are chain or 
hyphen galaxies. Hardly any DDGs were edge-on disks, 
which is understandable since these are most likely to 
have detectable rotation. Mergers are not a large frac- 
tion of either RDGs or DDGs. 

Restricting the sample to z > 0.4 changes the ratios 
little, other than to decrease the number of edge-on disk 
RDGs. There is a definite offset in morphology between 
RDGs and DDGs, but also substantial overlap. These 
morphologies are subjective; the link between kinematic 
and photometric morphologies will be addressed fur- 
ther with larger samples and objective measures in the 
DEEP2 survey. 

A significant number of galaxies were classed as chain 
or hyphen, types which are rare at low redshift. Two- 
thirds of the chain/hyphen galaxies were DDGs; only 
one-third of these elongated objects were dominated by 
rotation. TKRS 5627 shown in Figures and 01 is an 
example of a chain galaxy, elongated but nonrotating. 
Similar objects have been found at z ~ 2 by Erb et al. 
(2004): their slit spectra of 13 elongated galaxies showed 
dispersion, but rarely rotation. These results suggest 
that kinematic surveys must use caution in interpreting 
ellipticity: elongated objects may not be inclined disks 
and may not be rotating. 

The higher fraction of peculiar objects at high redshift 
suggests that galaxies evolve in photometric morphology, 
with a greater fraction of seemingly irregular galaxies at 
z ~ 1 (e.g. Abraham et al. 1996; Lotz et al. 2006). 
Possibly they also evolve in kinematic morphology. Lu- 



minous dispersion dominated galaxies appear at z ~ 1 
but are uncommon locally; even faint local blue galax- 
ies are usually dominated by rotation (e.g. Swaters et 
al. 2002). Perhaps DDGs settle from a disordered or 
non-planar kinematic state into more orderly rotation; 
or their disordered kinematics are dominated by bright 
star forming regions which then fade to allow a more or- 
dered background to be seen; or they evolve in luminosity 
by fading more than the overall population; or they cease 
forming stars and both fade and drop out of emission-line 
kinematic samples, leaving behind dispersion-supported 
faint red galaxies. 

7. CONCLUSIONS 

We have measured line-of-sight kinematic linewidths 
a\d for ~ 1000 galaxies in the Team Keck survey of the 
GOODS-N field and spatially resolved line-of-sight rota- 
tion (Vrot) and dispersion {a 2d) profiles for 380 of these. 
Most galaxies with linewidths are on the blue side of the 
color bimodality. Linewidths from integrated spectra are 
a measure of internal kinematics that is relatively ro- 
bust against observational effects, based on simulations 
of velocity fields and on comparisons of the linewidths 
with the spatially resolved kinematics. However, the un- 
known geometries of high-redshift velocity fields, scat- 
ter in galaxy ellipticities, and the lack of correction for 
projection into radial velocities mean that there is a fair 
amount of scatter between an individual galaxy linewidth 
measurement and a true circular velocity or dynamical 
mass estimate. For rotating galaxies, a±d ~ 0.6V C in 
the mean, due to a combination of the inclination and 
the fact that aid is less than line-of-sight V ro t\ this fac- 
tor should be kept in mind when estimating dynamical 
masses. 

The rotation and dispersion profiles of spatially re- 
solved galaxies show that not all galaxies exhibit a 
conventional rotation curve. Galaxies can be roughly 
divided into rotation and dispersion dominated, with 
V ro t > &2<i or vice versa. Dispersion dominated galaxies 
are blue, mostly irregular, and are not ellipticals; they 
probably have an effective dispersion, from disordered 
kinematics which are integrated over by the seeing. Dis- 
persion dominated galaxies exist at all redshifts in our 
sample, but low-z DDGs are quite faint. 

When line-of-sight rotation and dispersion are com- 
bined to make an estimate of the overall dispersion of 
the potential, as in S0.5 2 = 0.5V rot 2 + a 2 d 2 , the result 
correlates well with the integrated linewidth aid, demon- 
strating both that aid is a robust velocity indicator and 
that it is possible to construct scaling relations with ve- 
locity for a population of diverse kinematic properties. 
In Paper II we use the linewidths to measure evolution 
in the Tully-Fisher relation. 

At z ~ 1, rotation and dispersion dominated galaxies 
have similar magnitudes and colors, but somewhat dif- 
ferent size and morphology distributions. Rotation dom- 
inated galaxies are mostly disk and irregular types; dis- 
persion dominated galaxies are in the mean smaller in 
physical half-light radius, and ~ 60% of DDGs are ir- 
regular and chain/hyphen types. Only a small fraction 
of either RDGs or DDGs are obvious mergers. About 
two-thirds of chain and hyphen galaxies are dispersion 
dominated: elongated high-redshift objects cannot be as- 
sumed to be inclined rotating disks. It is not clear what 
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the local counterparts of dispersion dominated galax- 
ies might be; integral field spectroscopy, especially with 
adaptive optics, may shed some light on their nature at 
high redshift. 

There is some evidence that dispersion dominated 
galaxies are rare at bright magnitudes at low redshift, 
just as previous studies have shown that irregular pho- 
tometric morphologies in bright galaxies are more com- 
mon at high redshift than locally (e.g. Abraham et 
al. 1996; Lotz et al. 2006). It is possible that galax- 
ies evolve in kinematic morphology, settling into ordered 
rotation, or that luminosity evolution causes them to ap- 
pear more ordered, or to drop out of bright emission-line 
samples. These and other possibilities are avenues for 
further study with larger samples and quantitative mor- 
phological measurements. 
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